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Q ■ Abstract 

o 

We investigate the covariant Hamiltonian symplectic structure of General 

"o '. 

Q^i Relativity for spatially bounded regions of spacetime with a fixed time-flow 

bX). vector field. For existence of a well-defined Hamiltonian variational principle 



taking into account a spatial boundary, it is necessary to modify the standard 
Arnowitt-Deser-Misner Hamiltonian by adding a boundary term whose form 
depends on the spatial boundary conditions for the gravitational field. The 
most general mathematically allowed boundary conditions and corresponding 
boundary terms are shown to be determined by solving a certain equation 
obtained from the symplectic current pulled back to the hypersurface bound- 
ary of the spacetime region. A main result is that we obtain a covariant 
derivation of Dirichlet, Neumann, and mixed type boundary conditions on 
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the gravitational field at a fixed boundary hypersurface, together with the 
associated Hamiltonian boundary terms. As well, we establish uniqueness of 
these boundary conditions under certain assumptions motivated by the form 
of the symplectic current. Our analysis uses a Noether charge method which 
extends and unifies several results developed in recent literature for General 
Relativity. As an illustration of the method, we apply it to the Maxwell 
field equations to derive allowed boundary conditions and boundary terms for 
existence of a well-defined Hamiltonian variational principle for an electro- 
magnetic field in a fixed spatially bounded region of Minkowski spacetime. 
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I. INTRODUCTION 



The mathematical structure of General Relativity as a Hamiltonian field theory is well- 
understood for asymptotically flat spacetimes. As first shown by Regge and Teitelboim 
[1], with asymptotic fall-off conditions on the metric, there is a modification of the stan- 
dard Arnowitt-Deser-Misner (ADM) Hamiltonian [2] whose field equations obtained from 
the Hamiltonian variational principle yield a 3+1 spht of the Einstein equations. The ADM 
Hamiltonian itself yields the 3-1-1 Einstein equations only if compact support variations of 
the metric are used in the variational principle. For metric variations satisfying asymptotic 
fall-off conditions, the ADM Hamiltonian does not give a well-defined variational princi- 
ple since its variation produces asymptotic boundary terms that do not vanish. However, 
the boundary terms can be canceled by the addition of a surface integral term at spatial 
infinity to the ADM Hamiltonian. The resulting Regge-Teitelboim Hamiltonian yields a 
well-defined variational principle for the Einstein equations with asymptotic fall-off condi- 
tions on the metric. On solutions of the Einstein equations the Hamiltonian reduces to a 
surface integral expression over spatial infinity, which turns out to yield the physically impor- 
tant ADM definition of total energy, momentum, and angular momentum for asymptotically 
fiat spacetimes. Various modern, covariant formulations of this Hamiltonian structure are 
known [3-7]. 

A natural question to investigate is whether this Hamiltonian structure can be extended 
to spatially bounded regions of spacetime. An important motivation is astrophysical appli- 
cations where asymptotically fiat boundary conditions are not appropriate, e.g. collapse to 
a black-hole, mergers of binary stars, or collision of black-holes. Another important appli- 
cation is for numerical solution methods of the Einstein equations. In these situations the 
spatial boundary is not an actual physical boundary in spacetime, but rather is viewed as a 
mathematically defined timelike hypersurface whose boundary conditions effectively replace 
the dynamics of the gravitational field in the exterior region. 

In this paper and a sequel [8] , we work out the covariant Hamiltonian structure of Gen- 
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eral Relativity for arbitrary spatially compact regions of spacetime E x R whose spacelike 
slices possess a closed two-surface boundary 5E, with a fixed time-flow vector fleld tangent 
to the timelike hypersurface boundary x R. Rather than start with given boundary con- 
ditions on the metric, we instead seek to determine both the most general surface integral 
term necessary to be added to the ADM Hamiltonian in covariant form together with the 
most general corresponding boundary conditions on the metric at 9E such that the modi- 
fied Hamiltonian has well-defined variational derivatives. This would yield the most general 
mathematically allowed variational principle for the Einstein equations with spatial bound- 
ary conditions on the metric. To carry out the analysis we employ the covariant Hamiltonian 
formalism (referred to as the Noether charge method) developed in Ref. [4,9,10]. 

The main results are that we find Dirichlet and Neumann type boundary conditions for 
the metric at a spatial boundary two-surface and obtain the associated Hamiltonian surface 
integrals. Under some natural assumptions motivated by the symplectic structure arising 
from the ADM Hamiltonian, the most general allowed boundary conditions are shown to 
be certain types of mixtures of the Dirichlet and Neumann ones. We also investigate the 
geometrical structure of the Dirichlet and Neumann covariant Hamiltonians. These each 
turn out to involve an underlying "energy-momentum" vector at each point in the tangent 
space of the spacetime at the two-surface. In the Dirichlet case, this vector depends only 
on the extrinsic geometry of the spatial boundary two-surface. Most strikingly, when the 
vector is decomposed into tangential and normal parts with respect to the two-surface, the 
normal part yields a direction in which the two-surface has zero expansion in the spacetime. 

In Sec. II we first apply the Noether charge formalism to investigate, as an illustrative 
example, the covariant Hamiltonian structure of the free Maxwell equations on spatially 
compact regions of Minkowski spacetime. We show that this analysis leads to Dirichlet and 
Neumann type boundary conditions on the electromagnetic field, corresponding to conductor 
and insulator type boundaries, as well as mixed type boundary conditions which are linear 
combinations of the Dirichlet and Neumann ones. We also investigate more general bound- 
ary conditions which give rise to a well-defined Hamiltonian variational principle for the 
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Maxwell equations, and we obtain a uniqueness result for the mixed and pure type Dirich- 
let and Neumann boundary conditions under some assumptions. The associated Dirichlet 
and Neumann Hamiltonians are shown to reduce to expressions for the total energy of the 
electromagnetic field, including contributions from surface electric charge and current due 
to the boundary conditions. 

In Sec. Ill we carry out the corresponding analysis of the covariant Hamiltonian structure 
of General Relativity for arbitrary spatially compact regions of spacetime with a closed two- 
surface boundary, without matter fields. We make some concluding remarks in Sec. IV. In 
an Appendix we develop the Noether charge method for general Lagrangian field theories 
with a time symmetry. This approach extends and unifies some aspects of the covariant 
Hamiltonian formalism introduced in recent literature [4,9-11]. (Throughout we use the 
notation and conventions of Ref. [12].) 

Inclusion of matter fields and analysis of the geometrical properties of the resulting 
Dirichlet and Neumann Hamiltonians for General Relativity will be investigated in Ref. 
[8]. It is also left to that paper to discuss the relation between these Hamiltonians and the 
Regge-Teitelboim Hamiltonian in the case when the two-surface boundary is taken in a limit 
to be spatial infinity in an asymptotically flat spacetime. 



II. ELECTRODYNAMICS 

To illustrate our basic approach and the covariant Hamiltonian (Noether charge) for- 
malism, we consider the free Maxwell field theory in 4-dimensional Minkowski spacetime 

We use the standard electromagnetic field Lagrangian, where the field variable is the 
electromagnetic potential 1-form A^, with the field strength 2-form defined as F^^, = d^^A^^y 
The Lagrangian 4-form for the field A^ is given by 

^abcdi^) = 2^abcd^mn^"''' = ^^[ab*^cd] (2-1) 
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where *F^j, = e^f, is the dual field strength 2-form defined using the volume form e^^^. 

A variation of this Lagrangian gives 

SKUA) = 5[a0M(^' + (2-2) 

where 

^UASA) = 6SA^,*F^^ (2.3) 
defines the symplectic potential 3-form. Prom Eq. (2.2), one obtains the field equations 

^bcM) = 6a[,*F,,] = 0, (2.4) 
or equivalently, after contraction with the volume form, 

d^F^b = d'^W] = (2.5) 

which is the source- free Maxwell equations for A^. 

Let = [d/dtY be a timelike isometry of the Minkowski metric, with unit normalization 
C^Tjgjj = — 1, and let E be a region contained in a spacelike hyperplane t — Q orthogonal to 

with the boundary of the region being a closed 2-surface 9E. Denote the unit outward 
spacelike normal to 9S in S by s", and denote the metric and volume form on by 
<^a6 = ^ab - ^ah + ^ud = ^abcd^''^'^' Let and 9Et be the images of E and 9E 
under the one-parameter diffeomorphism generated by on Minkowski spacetime. Denote 
the metric compatible derivative operator on 9E by T)^. Let 

A^a = CdAa + K^aC = '^Cd^Aa^ + ^a(r^e)> (2-6) 

which is the Lie derivative of A^ with respect to 

The Noether current 3-form associated to ^" is given by 

JabM^ A) = ©a6c(^. A^) + "^^'L^bcai^) = 6*F[,^£,^„, + 12e'^F[„,*F^^, (2.7) 
which simplifies to 
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JaJ^, A)- 69[J*F,,je'^AJ + 2r6„,^(5fF_F-" - AF^^f'^) - CA^S^JA) (2.8) 

after use of Eqs. (2.4) and (2.6). Hence, one obtains the Noether current on solutions A^, 

Ja,M^ A) - (id^ai*F,ci^%) + 2rw(<5fF_F-" - 4F,„F^"). (2.9) 

(Note, one easily sees that this 3-form Jabd^^A) is closed but is not exact, i.e. there 
does not exist a Noether current potential Q^^^^A) satisfying Jabd^^A) — 3d^g^Q^^^{^, A).) 
Correspondingly, the Noether charge on solutions A^ is given by 

QdO - X ^^"^^^-^ ^) = L ^«^e/^(-8^en^'" + ^^tF^nF^ + £ 2*F,/A,. (2.10) 

This expression simplifies in terms of the electromagnetic stress-energy tensor defined by 

Tf(F) = 2F,„F^- - \stF^^F-. (2.11) 

Thus, 

\Qi:iO= j^Ci'T^{F)d'x + j^^CA^es'F^dS (2.12) 

where d^x and dS denote the coordinate volume elements on E and obtained from the 
volume forms Cabcd^'' ^&c' respectively. 

A. Covariant Hamiltonian formulation 

The symplectic current, defined by the antisymmetrized variation of O^cdl^' ^A), is given 
by the 3-form 

\^bcdi^iA, S^A) = 5^A^,52*F^ - 52A^,5,*F^. (2.13) 

Then the presymplectic form on E is given by 

n^{6^A,62A) = J^uj^^{6,A,62A). (2.14) 

A Hamiltonian conjugate to ^ on E is a function Hy,(^; A) = /s^abc(^' A) for some locally 
constructed 3-form H„^p(^; A) such that 
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5H^(^; A) = H^i^; A, 5A) = n^{5A, C^A) (2.15) 

for arbitrary variations 5A^ away from solutions A^. 

Prom the expression (2.7) for the Noether current, the presymplectic form yields 

n^{SA,C^A)= JjJa,M,A)+At5A^,S^,^^{A) - £re„,,(A, 5A). (2.16) 

Hence, for compact support variations 5A^ away from solutions A^^, the Noether current 
gives a Hamiltonian (2.15) with 'H^(,^(^; A) = J^j^^iC, A), up to an inessential boundary term. 
The simphfied expression (2.8) for J^bd^^ ^) thereby yields the Hamiltonian 

^(C; ^) = 4 / ei\T^{F) + A^d^Fjd'x. (2.17) 

On solutions A^, this Hamiltonian is equal to the total electromagnetic field energy on E, 

To define a Hamiltonian (2.15) for variations 5A^ without compact support, it follows 
that the term C^abd-^^ ^■^) (2-16) needs to be a total variation at the boundary 9E, 

i.e. there must exist a locally constructed 3-form B^^J^A) such that one has 

f e>a6c(^, <^^)laE = {C^B^^iA) + d^^a,^{i- A, 5A))U (2-18) 

where af^{^]A,5A) is a locally constructed 1-form. This equation holds if and only if, by 

taking an antisymmetrized variation [10], one has 

e^e^aboi^i^^ 52A)|as = 5,A, 52A)\aT. (2.19) 

where 

d^A) = e''6,a,{i- A, d^A) - e'^'d^a.^i; A, 6,A) (2.20) 

is a locally constructed vector, in T(9E), which is skew bilinear in 5iA,62A. The term 
involving the symplectic current is given by 

'Tu^a^ci^^A, <^2A)= 9>s%''\5,A^5,F^, - S^AJ^FJ (2.21) 
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with /i^j^ = rj^j^ — s^Sj, = (7^^ — Given a solution of equation (2.19), one can then 

determine B^f^^{A) from equation (2.18) by 

= Ss%'^FJA^ - V^a^iC; A, SA) (2.22) 

where 

a''{^;A,5A) = e'''a,{^;A6A). (2.23) 

This leads to the following main result. 

Proposition 2.1. A Hamiltonian conjugate to on E exists for variations SA^ with 
support on if and only if 

8h'"'s'{S^A,^^J,A^^ - 52A,^^J^A^^)\^^ = Vj'^iC; S,A, S,A) (2.24) 

for some locally constructed vector l3g^{^,diA,d2A), in T[dYl), which is skew bilinear in 
5iA,52A. The solutions of equation (2.24) of the form 5T^{A)\q^ — give the allowed 
boundary conditions J-'^{A)\qy: for a Hamiltonian formulation of the Maxwell equations in 
the local spacetime region E^, t > 0. For each boundary condition, there is a corresponding 
Hamiltonian given by the Noether charge plus a boundary term 

H^iC; A) = / J„,,(e, A) - I tB^UA) ^ H{C; A) + Hs{C; A), (2.25) 

with 

H{i- A)=aJ^ ei\T^{F) + A^d'^Fjd'x, (2.26) 

H^{i- A) = £ Ci^A^i's^F,^ \B^^{A))dS, (2.27) 

where B^[A) — e"^B^^J^A) is determined from 

iCBM) - Ss%''FJAJ\ej: = V^a\C;A,5A)\ej: (2.28) 

with a'^{^] A,6A) given by Eqs. (2.20) and (2.23). Note, B^[A) is unique up to addition of 
an arbitrary covector function of the fixed boundary data J^^(A). 
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The results in Proposition 2.1 take a more familiar form when expressed in terms of the 
electric and magnetic fields on E defined by = 2^^F^^, = C^*-^^^, which are vectors in 
T(E) (i.e. ^"'Eg^ = CBg^ = 0). A convenient notation now is to write vectors in T(E) using 
an over script — >, and for tensors on M, to denote tangential and normal components with 
respect to E by subscripts || and ±, and denote components orthogonal to E^ by a subscript 
0. Then we have 

E = dAQ- d^A, B = dxA. (2.29) 

In this notation, the presymplectic form (2.14) and the Hamiltonian (2.25) reduce to the 
expressions 

\nj:{5iA,52A)^ [ {52A-5iE-5iA-62E)d^x 

= J^{6iA ■ {d^^S^A - dS2A^) - S2A ■ {d^S.A - d6iA^))d^x (2.30) 

and 

^-H^it, A) = ^{E' + B') + A,d ■ E^ d'x - £ A,E^ + \B^{A)dS. (2.31) 

Note that the Hamiltonian field equations obtained from Hs{^;A) are given by the 
variational principle 

H^{^; A, 6 A) = nj:{6A, C^A) (2.32) 

for arbitrary variations ^^als- These field equations split into dynamical equations and con- 
straint equations, corresponding to a decomposition of A^ into dynamical and non-dynamical 

— * 

components, respectively A and Aq, determined by [4] the degeneracy of the presymplectic 
form (2.30). In particular, this yields the Gauss-law constraint equation 

d-E = AAQ-dod-A = (2.33) 

obtained from H^{^] A,6Aq) = through variation of Aq, and the dynamical Maxwell 
evolution equation 
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d^E-dxB^ {-d^^ + A)A + d{d^^A^ - 9 • A) = (2.34) 

— * — * — * — * — * 

obtained from H!^{^\A^5A) — Q,y,\5A^C^A) through variation of A^ where A = 9 • 9 is 
the Laplacian on E. Thus, the Noether charge (covariant Hamiltonian) formahsm here is 
equivalent to the standard canonical formulation [12] of the Maxwell equations. 

B. Dirichlet and Neumann Boundary Conditions 

Two immediate solutions of the determining equation (2.24) with = are boundary 
conditions associated with fixing components of A^ or F^^ = d^^^A^-^ at dJ^t for t > 0. 
Consider 

(D) a.'^MJas, -0, f^AJas, = 0, t>0, (2.35) 
or equivalently 6A\\ = 6A^^ = 0, for t > 0, called Dirichlet boundary conditions, i.e. 

•^a (^) = (2.36) 

(N) <7,V5F„JaE, = 0, s^eSFja^^^O, t>0, (2.37) 

or equivalently dj_ x SA\\ — d^^ x SA± = dj_5AQ — dQdA± = 0, for t > 0, called Neumann 
boundary conditions, i.e. 

•^a (^) = K'^'F,^- (2-38) 

Theorem 2.2. For the boundary conditions (D) or (N), a Hamiltonian conjugate to ^" 
on evaluated on solutions A^, is given by 

l//°(e; A)^ 2 e.^.M'F^nF''' - ^^'F^nFn + 2 £ e^^CA^^'s^F,^ (2.39) 

= ^ / E'^ + B'^ d^x - i AqE^ dS, (2.40) 

li/^(e; A)= 2 e.^U'^CF^^F'- - I^'F^F^) + 2 £ e,,a^^A,/F,, (2.41) 

= ^ / E^ + B^ d^x- I iA X B)x (2.42) 

2 7s JdT, 
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Proof: 

For (D), one has 



= 8(f - ^'"'SAyrj = (2.43) 



and hence = 0, so thus CB^{A) = 0. For (N), one has 



e'rQai>M,jC,A)^6e'r{SA^a*Fbc]) 



S{6e'rA^^*F,^^) + Si-CA^^s^SF,^ + A^VSEJ (2.44) 



with B^,^{A) = 6A[„*F,,j. Thus, 

CBM) = 66^r^[,*F,,j = 8s,aX\ (2.45) 



□ 

Note that for both boundary conditions (D) and (N), the surface integral terms in the 
Hamiltonian take the form 

HB{i;A)^Ai ePSA)dS (2.46) 

where 

A^es'F,^, (2.47) 
PliA)= -i/a'^A,F^. (2.48) 

There is simple physical interpretation of the (D) and (N) boundary conditions: (D) 

— * 

involves fixing A^ and Ay at dT^t for i > 0, which means 

^11 = - dJw, = 9,1 X A\\ (2.49) 
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are specified data at the boundary surface (analogous to a conductor) as a function of time. 

— * — * 

Note, consequently, E± and S|| are left free by the boundary condition (D) and therefore are 
induced data for solutions Aq,A of the Hamiltonian field equations; (N) reverses the role of 

— * — * 

the induced and fixed data at the boundary surface, so now E± and B\\ are specified as a 
function of time (analogous to an insulator), while the induced data for solutions A^^, Aof the 
Hamiltonian field equations are E\\ and B±. Note the fixed data here are gauge-equivalent 

— * 

to specifying the normal derivative of Aq and at dUf in Et for i > 0, 

E^^d^{A,-d,x), S||=9^x(l||-9||X), (2.50) 

— * — * 

where x is given by dj^x — 

Moreover, the Hamiltonians (2.40) and (2.42) on solutions Ag,Ah.eiye the interpretation 
as expressions for the total energy of the electromagnetic fields, with the surface integral 
parts representing the energy contribution [13] from an effective (fictitious) surface charge 
density in the case (D), and effective (fictitious) surface current density in the case (N), 
associated to the specified data at the boundary surface. In particular, effective surface 

— * — * 

charges and currents arise, respectively, when E± or Ay are left free [13] on the boundary 
surface. 

We remark that similar boundary terms arise in the asymptotic case when the boundary 
surface is taken to be a 2-sphere at spatial infinity on E (see Ref. [14]). 

C. Determination of allowed boundary conditions 

The symplectic current component e"^^"'ujg^i^^{diA, 52A) involves only the field variations 
hJ'6A^, hJ's'^d^^SA^y We refer to the components o-^^A^,^'^A^,a^'^s°'F^^,^'^s°'F^^, or equiva- 
lently A^, E±, as symplectic boundary data at 9E. Hence, in solving the determining 
equation (2.24) for the allowed boundary conditions on A^, it is then natural to restrict 
attention to boundary conditions involving only this data. (Some remarks on more general 
boundary conditions are made at the end of this section.) To proceed, we suppose that the 
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possible boundary conditions are linear, homogeneous functions of the symplectic boundary 
data, with coefficients locally constructed out of the geometrical quantities s^, a^^, e^^ 
at the boundary surface. We call this type of boundary condition a symplectic boundary 
condition. 

Theorem 2.3. The most general allowed symplectic boundary conditions 

^.iK'^a, hfs'F,^; f , /, ej (2.51) 
for existence of a Hamiltonian conjugate to ^" on E are given by 

= bo'7,'A^ + ao<«"F„e + bi^bCA^ + ^M'^^^^ca^ (2-52) 

or equivalently, 

a^iboSA^ + ao/(5F„J |as* = 0, t> 0, (2.53) 

C{biSA^ + ais'^^F J las, =0, t > (2.54) 

for any constants ao, 6o (not both zero), oi, bi (not both zero). 
Proof: 

First, we show that = without loss of generality in the determining equation (2.24) 
for boundary conditions of the form (2.51). Note the left side of Eq. (2.24) is algebraic in 

5F^^ = d^^SA^^y Since the right side necessarily involves at least one derivative on 6A^, 
we must have 

= 6,A,62A^(3''"' (2.55) 

for some tensor 

^abc ^^a^elbc]_ ^2.56) 

We substitute expression (2.55) into Eq. (2.24) and collect all terms that do not involve 

just the symplectic boundary data, namely s^A^, C^^d^aci ^b^d^aci ^^"^ ^(a^b)- '^^^ 
coefficients of these terms yield algebraic equations 
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g^^H)c^O^ g^^Mc^O, = 0. (2.57) 

Then, since (3"'^'^ has the form (2.56), we find that the solution of the equations in (2.57) is 

= (2.58) 

and so = 0. 

Hence, the determining equation (2.24) reduces to 

h!"'s\5,A,d^^5,A^^ - 5i A9[a<^i^c]) laE = 0, (2.59) 

which we are now free to solve as a purely algebraic equation in terms of the variables 
and 5F^ = d^JA^^,\.e. 

h%\S,A,S,F^^ - 5,A,5,FJ = 0. (2.60) 

It is straightforward to show from the form of Eq. (2.60) that the only solution which is 
linear, homogeneous in the previous variables is given by 

Hi.^M, = n2,^/5F„, (2.61) 

where IIi'"^ and 112'"^ are some symmetric tensors orthogonal to s", such that Eq. (2.61) can 
be solved for either ^^SAj^ or ^^s"'dF^i^, and either a^5A^ or cr^s'^SF^i^. Since we require the 
coefficients in the boundary conditions under consideration to be locally constructed out of 

C", s", cTftc' ^ftc we see that 

Ui'^^boa' + hC'e, U2"'^aoa'' + a,C'e, (2-62) 

for some constants oq, oi, bo, hi with oq 7^ or &o 7^ 0, and oi 7^ or 61 7^ 0. This yields the 
general solution (2.53) and (2.54) given in the Theorem. □ 

The boundary conditions given by Theorem 2.3 comprise the following separate types: 
(i) for ao = ai = or 60 = ^'i = 0, one obtains, respectively, Dirichlet (2.35) and Neumann 
(2.37) boundary conditions; (ii) for 60 = 0(ao 7^ 0),6i 7^ 0, the boundary conditions yield a 
one-parameter ai/61 = ci family of the form 
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(JfeV^F^JaE, = 0, fMJaE, = -cisT<^F„Jas„ ^ > 0, (2.63) 

or equivalently, 5B\\ — s5Af^+ci^SE± — 0, for t > 0; (iii) similarly, for bi — 0(ai 7^ 0), 60 7^ 0, 
the boundary conditions yield another one-parameter ao/bo = Cq family of the form 

<<^Aj5s. = -cos"5i=^„,|aE„ es''SFje^, = 0, t>0, (2.64) 

or equivalently, s x 5A\i + co^6B\\ = 5E± — 0, for t > 0; (iv) finally, for bo 7^ 0, 61 7^ 0(ao 7^ 
0, ai 7^ 0), we obtain a two-parameter ai/61 = ci and ao/60 = cq family of the boundary 
conditions 

(^bSAjQj:^ = -CqS^'SFJqj:^, C^A^ldSt = -Cis''CSFjoj:^, t > 0, (2.65) 

or equivalently, sx 5A\^ + co^5B\\ — s5Aq + ci^5E± — 0, for t >0. 

The fixed data for these boundary conditions (2.63) to (2.65) corresponds to specifying, 
respectively, the field components 

sA^ + C^E^, B\\, CiT^O, (2.66) 

E^, sx A||+co^S||, Co 7^0, (2.67) 

1 - - 1 - 

sA^ + Ci-E^, sx^ll + co-5||, co7^0,ci7^0, (2.68) 

at the boundary surface for t > 0. From Proposition 2.1, we readily obtain the Hamil- 
tonian boundary terms corresponding to these boundary conditions (by a proof similar to 
that for Theorem 2.2). 

Theorem 2.4. For boundary conditions (2.63) to (2.65), there is a respective Hamiltonian 
(2.25) conjugate to ^" on with boundary terms given by 

HbH; A)^ 4 £ e^^iCA^es^F,^ + a^'A^F,^ - c^i^s^Fj^) 

= - / 2(1 X B)s_ + ci^x^ dS, (2.69) 

Jdll 

HB{i;A)^^j^^e,^ca'^s^F^s-F^^ 

coB\^ dS, (2.70) 



as 
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= - 2A^Es_ + ciE^^ - co^||2 dS. (2.71) 

Interestingly, among the allowed boundary conditions given by Theorem 2.3, we observe 
that the Hamiltonian boundary terms vanish identically in one (and only one) case, when 
Co = in boundary condition (2.64), i.e. 

<5^JeE, = 0, eV5F„>s, = 0, t>0, (2.72) 

or equivalently, dA\^ = SE± = 0, for t > 0. The resulting Hamiltonian (2.25) reduces, on 
solutions A^, simply to the expression for the total energy of the electromagnetic fields, 
H{^;A) = 4Jj.CtTae{F)d^x = fj^E^ + (fx. The fixed data corresponding to the 
boundary condition (2.72) is A\\ and E±, which means that the normal components of the 
electric and magnetic fields at are specified for i > 0, 

= 9|| X A\\, E± = d^A^ - doA±. (2.73) 

— * — * 

Since A\\ and E± are fixed, there are no effective charges and currents associated to the 
boundary surface. Thus, the total electromagnetic energy involves no surface integral con- 
tributions in this case. 

D. Remarks 

We conclude with some short remarks on uniqueness of the boundary conditions obtained 
in Theorem 2.3. 

Note that the symplectic boundary conditions (2.53) and (2.54) are hnear combinations 
of the tangential and normal parts of the Dirichlet and Neumann boundary conditions, 
referred to as mixed boundary conditions. In physical terms, they correspond to specifying 
boai^^A^ + aos"'af^'^F^^ and bi^'^A^ + ai^'^s"'F^^ as boundary data at dT,t for t > 0. Theorem 2.3 
gives a uniqueness result for these mixed boundary conditions under the natural assumption 
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(2.51) about the general type of boundary condition considered on the fields at the boundary 
surface. If this assumption is loosened, then there exist additional boundary conditions 
allowed by the determining equation (2.24). 

In particular, one can trade off some of the mixed boundary conditions on the sym- 
plectic boundary data for boundary conditions involving the symmetrized derivatives of 

at 9E. For example, an allowed boundary condition satisfying Eq. (2.24) is given by 
^(A) = cr^'^d^Af,, a^-^s^d^f^A^^), or equivalently 

eSAJe^,^0, = 0, a^'s%SA^^\9j:, ^ 0, t>0, (2.74) 

with = Se^^A^s^SA^. Prom Eq. (2.28) one obtains 

ICB^ = s'a'\A^d^Ae) - KdAa) - A^'«^ (2-75) 

where T>^s^ — crj^crj^d^s^ is the extrinsic curvature of in E. Hence the corresponding 
boundary term in the Hamiltonian is given by Eq. (2.46) with 

Pa(A) = A^'CF,, - s'a'^iA^d^.A^^ ^ A^dA,) + ^A.A^P'^/. (2.76) 

III. ANALYSIS OF GENERAL RELATIVITY 

We now apply the Noether charge analysis to General Relativity, specifically to the 
vacuum Einstein equations for the gravitational field in a spatially bounded spacetime region 
with a fixed time-fiow vector field. It is straightforward to also include matter fields, as we 
discuss in Ref. [8]. 

For General Relativity without matter sources, the starting point is the standard La- 
grangian formulation of the Einstein equations with the spacetime metric as the field vari- 
able. It turns out, however, that the analysis is considerably simplified by introduction of 
a tetrad (orthonormal frame). Moreover, taking into account local rotations and boosts 
of the tetrad, the boundary conditions and resulting Hamiltonians that arise in the tetrad 
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formulation are equivalent to those obtained purely using the metric formulation, up to a 
boundary term in the presymplectic form. 

After setting up some preliminary notation and results in Sec. Ill A, we will first consider 
a Dirichlet boundary condition as explained in Sec. IIIB. Then we will carry out details of 
the Noether charge analysis with the Dirichlet boundary condition using the tetrad formu- 
lation of General Relativity in Sec. Ill C. The resulting covariant Dirichlet Hamiltonian for 
General Relativity is summarized in Sec. Ill D where we will discuss the equivalence between 
the metric and tetrad formulations. In Sec. Ill E we will investigate a Neumann boundary 
condition and corresponding Hamiltonian, along with more general boundary conditions and 
Hamiltonians. The main result will be to establish uniqueness of mixed Dirichlet-Neumann 
boundary conditions for existence of a Hamiltonian formulation of General Relativity. Fi- 
nally, in Sec. Ill F we will briefly discuss the form of the Dirichlet and Neumann covariant 
Hamiltonians, and relate these to an analysis of boundary terms for the ADM Hamiltonian 
using the standard (non-covariant) ADM canonical variables. 



On a given smooth orientable 4-dimensional spacetime manifold M, let g^^^ be the space- 
time metric tensor, ^ahcdid) ^® volume form normalized with respect to the metric, and 
^V^ be the covariant (torsion-free) derivative operator determined by the metric. 

Now, let be a complete, smooth timelike vector field on M , and let E be a region con- 
tained in a spacelike hypersurface with the boundary of the region being a closed orientable 
2-surface (9S. Let denote the unit outward spacehke normal to orthogonal to let 
denote the unit future timelike normal to dTi orthogonal to s"" . Denote the metric tensor 
and volume form on by 



A. Preliminaries 
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ab 



9ab 



— ss, -\- tt,, 
a a 0' 



(3.1) 
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abed 



(^)sV. 



(3.2) 



This yields the decompositions 
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9ab = '^ab + ^ah " *a^6' 



(3.3) 
(3.4) 



Note that one has 

= ATf + AT", N^s^ = Cs^ = (3.5) 

for some scalar function N and vector function N"" on 9E. It is convenient to extend the 
previous structures off as follows. Let V be the spacetime region foliated by the images 
of E under a one-parameter diffeomorphism on M generated by and let B be the timelike 
boundary of V foliated by the images of 9S. Fix a time function t which is constant on 
each of the spacelike slices diffeomorphic to E under ^" in V and which is normalized by 
C9j> — 1, such that t — Q corresponds to E. Then B is a timelike hypersurface in M whose 
intersection with spacelike hypersurfaces given by t = const in V consists of spacelike 
2-surfaces SEj diffeomorphic to 9E. Finally, let s", f' , a^^, e^^, N , N"' be extended to 9Et, 
and let denote the unit future timelike normal to E^ parallel to dj,. 

Note that, by construction, is hypersurface orthogonal to B and hence 

^[A^a] = 0- (3.6) 

If is expressed as a linear combination of s^, c?^t, then since dJ: obviously is hypersurface 
orthogonal to dT^, it follows that 

SyatAU = 0. (3.7) 

In addition, note that s^dj, measures the extent to which E^ fails to be orthogonal to B. 

Let {VdE)a and {Vt)^ be coordinate projection operators onto the respective tangent 
spaces of the 2-surface dT^t and the integral curve of t"", and let {VB)a = {'PdT,)a + i'Pt)a, 
which is the projection operator onto the tangent space of the timelike hypersurface B. 
Note that these operators are independent of the spacetime metric, as they involve only the 
manifold structure of B and dUt in local coordinates in M. 
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Hereafter we work in terms of an orthonormal frame 9'^ (i.e. tetrad) for g^j^ on M. The 
frame components of s", f, n", a^^, e^^, g^^, e„,,cd(f ) are given by 

= a"'^^^^, e'^'^ = e^'ey,, (3.9) 

where r]'^" = a'^" + s'^s" — t'^t" = diag{—l, 1, 1, 1) is the Minkowski frame-metric, with 
a'^'^ — diag{0, 0, 1, 1). This leads to an orthonormal frame for the metric a^, given by 

= (3-11) 

satisfying 

sa/ = t,a/ = Q. (3.12) 
Let the inverse orthonormal frame for g^^ and for cr^^ be denoted by 

r.a ah nil a ah a /oio\ 

Then, one has the decompositions 

0: = ^a' + sy-tX, 9; = a; + s^s^-t\. (3.14) 

For later use, we will partially fix the S0(3,l) local gauge freedom in 6*^ by choosing 
the coefficients s^, in the frame decomposition (3.14) to be fixed functions on M, so that 
under a variation Sg^^^^, 

S^, = ^t^ = 6a^,^0 (3.15) 

and hence, correspondingly, 

K = s,se:, K = t,K^ (3-16) 

^<^a, = '^%r^J()l) = 2a(/5<7,)^. (3.17) 
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Similarly, one then also has 

KUf^) = K,,.(^la&yjO''^ = ^aU&WM, (3-18) 

= 2v^yC] = (3.19) 

and thus 

= = 0. (3.20) 
Consequently, some useful identities are given by 

K = ^'^f + - t'K^ (3.21) 

K = -V.S^c = + - (3-22) 
Now, a variation of the spacetime metric Sg^^^ can be decomposed into the parts 

<^^a6 = K6 + 2«(aS-2M^^>)- (3.23) 
By hypersurface orthogonality, one has the identities 

5s^ = sJSs,, 5t^ = sj5t, - t/5t, (3.24) 

and 

^'^ab = '^a^b'^^'^cd + ^aH^'^^^cd " ^aX^'^^^cd' (3-25) 

Then, from the relation 

5g''' = Sa""' + 2s("5/^ - 21^^" 51^^ = -g'"'g'"'5g^, (3.26) 
it straightforwardly follows that 

(5s" = '^h^s' + s"'Sj^5s' — t"'tfjSs', (3.28) 
5f = V^i'' - t\5t\ (3.29) 
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where 

Sj^Ss^ — —s'^Ss^^, t^^Ss'' — — i^(5sjj, t^^St'' — —t^6t^^, (3.31) 
and again by hypersurface orthogonality, 

sja"'' = tja"'' = 0, sjf = 0. (3.32) 
Thus, the hnearly independent parts of S9'^, or equivalently of 5^^, are given by 

SaJ', (T^5s' , f^b^t^i ^b^^^ 1 ^b^^^ 1 h^t'^- (3.33) 
Throughout, the time-flow vector field is taken to be fixed, S^"" — 0, under variations of 

B. Dirichlet boundary condition 

There is a natural motivation for a Dirichlet boundary condition on the gravitational 
field in the Einstein equations in analogy with the Maxwell equations where the tangential 
components of the electromagnetic field potential are specified at the boundary. For 
General Relativity, similarly, one can introduce a Dirichlet boundary condition given by 
specifying the tangential components of the spacetime metric g^f^ at the 2-surfaces dJ^f 
This boundary condition is expressed equivalently by conditions on the variation of the 
metric tensor 

SaJan^^O, 5tj9E, = 0, t>0. (3.34) 
Geometrically, this means that the metric given by 

Kb = ^ab - tJb (3-35) 
on the timelike boundary hypersurface B is specified data, so it is held fixed under variations 

of 9ab^ 
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5h^^ = on B. (3.36) 

The geometrical form (3.36) of the Dirichlet boundary condition is often introduced when 
one considers an action principle for General Relativity on a spacetime manifold with a fixed 
global timelike boundary hypersurface [12,15-17]. We will see in the next section that this 
boundary condition in the form (3.34) emerges naturally from the Noether charge analysis 
for the existence of a Hamiltonian formulation of General Relativity for a spatially bounded 
local spacetime region. 

Note that, from the relations (3.27) to (3.31), one can decompose the Dirichlet boundary 
condition (3.34) into an intrinsic part 

= -/i"'rt/i^|9s, = 0, i > (3.37) 

and an extrinsic part 

h^^,5s'\d^^ = -s'5h^^\d^^ = 0, t > (3.38) 

with respect to the timelike hypersurface B. The intrinsic part corresponds to fixing just 
the metric Vsh^^^ restricted to the tangent space of B, where the projection Vb removes 
components of the hypersurface metric proportional to s^. Correspondingly, note that the 
volume form 

on this surface is also fixed, 5&oji,JJi) — Q onB, since 

Kbcih) - le^.{h)h-Sh^^ = -le^cih)h^Jh-. (3.40) 

The Dirichlet boundary condition has a simple formulation in terms of the orthonormal 
frame 9'^. It is convenient to introduce a frame for the metric /i^^ by 

hS = hX = - tj", (3.41) 

and inverse frame 
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h""^ = (3.42) 
Then the Dirichlet boundary condition (3.34) is equivalent to 

<^/i/|as, = 0, i>0, (3.43) 

with intrinsic part 

^h''^^, = 0, i > 0, (3.44) 

which is equivalent to Eq. (3.37). These equivalences are immediate consequences of the 
identities h^^ — hj^hj^^rj^^ and h"''^ — h^hj^r]^^ . From these identities, one also has 

^K, = (3-45) 

K,c{h) = hfSh/. (3.46) 

An additional useful identity is given by 

5h: = -(/5.„ + sjs'^) = -sXSs' (3.47) 

and therefore Vs^h^ — 0. 

Finally, note that the intrinsic part of the Dirichlet boundary condition on the frame 
decomposes into 

= 0, ^rias, = 0, t > 0. (3.48) 

The full, extrinsic Dirichlet boundary condition is necessary and sufficient for 5h^\aT.t = 
0,t>0. 

C. Noether charge analysis 

We consider the standard tetrad formulation of General Relativity, using an orthonormal 
frame 9'^ for gi^^ and a frame-connection 

rne) = Q'^^vf, = 2^^['^a[„^^; - e'^-e-e^^d^.eiy (3.49) 
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Here the expression in the second equahty is obtained from the relation 

= VS- (3.50) 
The curvature of this connection (3.49) is given by 

Kno) = 2a[„r,j'^^(^) + 2r[/'^(^)r„/(^) = R^UgW^'o"-, (3.51) 

related to the Riemann curvature tensor Rabcdid) 9ab- 

With 9'^ as the field variable, the Lagrangian 4-form for General Relativity (without 
matter sources) is given by 

LabcM = ^aUmO) = QO^J:,Rc^,M (3.52) 

where 

KnO) = Rcaa.iOy""' = 2d^J,^'^ - r[/[''(e)f,j/l (9) (3.53) 

in terms of f„^^(^) = r„"^(^)e a0fj,u- Then the variation of L^^JJ)) gives, after integration 
by parts and use of the connection equation (3.50), 

where 

QUO. se) = i2e;9:5T,^^M = k,j:o;st:''{9) (3.55) 

defines the symplectic potential 3-form. The field equations for 9'^, obtained from the coef- 
ficient of 59'^ in Eq. (3.54), are given by 

^L(^) = 12^r.4.],.(^) = ^hc^mR'^'iO) - y'R{9)) = 0. (3.56) 

Thus 9'^ satisfies R"''^\9) = 0, which is equivalent to the vacuum Einstein equations for the 
spacetime metric 

Rjg) = (3.57) 
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arising as the stationary points of the action functional S{g) — Im ^abcdi9)-^i9) u^^der com- 
pact support variations of g^^^. 

The Noether current associated to is given by the 3-form 

= 12^;CAr„],.(^) + 24e'^fe^^^„,],.(^) (3.58) 

with the first term obtained from the variation of the frame connection (3.49) after replace- 
ment of 69'^ by the Lie derivative C^O'^ = C^e^a + ^e^a^^ f^^^ ^^^^ Lie 
derivatives commute with exterior (skew) derivatives. We now simplify the first term in 
Eq. (3.58) as follows. First we express 

Hence we obtain 

i2^;CAf„)^,(e) = d^^{i2e',e:fr^^^{e)) - i2r^;c4]e,.(^) (3.60) 

through use of the identity (3.53). Next we combine the second terms in both Eqs. (3.58) 
and (3.60) to get 

-12r^;C^aie..(^) + 24eX^.^^„,],.(^) = Ce^Mi^StRig) - 8Rf{g)). (3.61) 
Thus, one obtains the Noether current 

^a.c(e; 0) = 3a[„g,,] (C; e) - r^,/L(^) (3.62) 

where 

Qbci^;0)=^tra,MVc (3-63) 

is the Noether current potential 2-form. 

On vacuum solutions 9'^, the Noether current reduces to an exact 3-form 

JabM-^^)-^aQbc]i^-,&)- (3-64) 
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Therefore, the Noether charge for vacuum solutions is given by the boundary 2-surface 
integral 

QdO - 1 Ja^; 0)- H^^tr^.^ioy = ^^r.^Mt's'^ds (3.65) 

where dS is the volume element on corresponding to the volume form e^^ in local coor- 
dinates. 

Now the symplectic current, defined by the antisymmetrized variation of 0ftcd(^'^^)' 
given by the 3-form 

^^U^, 6,9) = 9^M,6,t,^^^i9) - 9^^^,9''J^f,^^^i9). (3.66) 

Then the presymplectic form on S is defined by 

n^{9,5AS29) = j^u:,J9,5,9,5,9) = 24 jj^M^.^^.i'^) - 9;5,9'',5,f,^^^{9). (3.67) 

A Hamiltonian conjugate to ^ on E is a function H-s{^;9) — /s'^a&c(^'^) some locally 
constructed 3-form 'Hg^^{C',9) such that 

5Hj:{C; 9) = H^i^; 9, 59) = n^{59, C^9) (3.68) 

holds for arbitrary variations 59^ away from vacuum solutions 9^. 

In terms of the Noether current (3.62), the presymplectic form on E yields 

Qs(e, 59, [ 5J^JC; 9) - 4t£^^^{9)59,^^ f ee,j9, 59). (3.69) 

Consequently, for variations 59'^ with compact support on the interior of E, the Noether 
current defines a Hamiltonian conjugate to 

which is equal to the Noether charge (3.65) when 9'^ is a vacuum solution. Explicitly, from 
Eqs. (3.56) and (3.63), one has 

0)- 8 e^aeO'iR;{9) - Ic'Rm + 4 e,^C%,Me"' 
= 8 Cn.iRm - \9:R{9))d'L + 8 j^^ ^Ta^Mt's'^dS (3.71) 
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where dT, is the coordinate volume element on E obtained from the volume form e^^^^n . 

To define a Hamiltonian Hj]{^; 9) for variations dO'^ without compact support, it follows 
that the term C^abd^^ (3-69) needs to be a total variation at the boundary 812, 

i.e. there must exist a locally constructed 3-form B^^J^O) such that one has 

re>a6c(^. = (f <^^„,c(^) + 5[a«;.](e; 0. m)\e^ (3.72) 

where q;(,(^; 9, 59) is a locally constructed 1-form. This equation is equivalent to 

^"'re>a6c(^. = e"'r<^^a.c(^. <^^)laE + 'y'd,a\^- 9, 59)\3^ (3.73) 

where a'^{^]9,59) — e^^ a^{^] 9 , 59) and the symplectic potential term is given by, using 
identity (3.4), 

e"'r0a.e(^, m= ^'^Cty^S^fcOl^^:''{0). (3.74) 

Hence we now have the following result. 

Proposition 3.1. A Hamiltonian conjugate to ^" on E exists for variations 59'^ with 
support on 9E if and only if 

e'^'eSB^.Ma^ = m^aS^fX^K^'imdi: - ^'d,a\i- 9, 59)\a^ (3.75) 

for some locally constructed 3-form B^^^{9) in T*(S) and locally constructed vector 
dt{^]9,59) in T((?E). The Hamiltonian is given by the Noether charge plus an additional 
boundary term 

M^;0)= [ja,Mo)-[eB^do) 

= 8 X Va^c(^) + ^ Qa,i^-, 0) - eB^,M 

= 8 Cn^iRfiO) - \0'em)d^ + C'i^^a,AO)t''s'' - lB,{9))dS (3.76) 
withB,{9)^e'%j9). 
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We now show that the equation (3.75) for existence of a Hamiltonian (3.76) is satisfied 
for the intrinsic Dirichlet boundary condition on 0^, 

^h''^^, = 0, t > 0, (3.77) 

and then we derive the corresponding Hamiltonian boundary term. Henceforth we take 
9^ to satisfy the gauge conditions (3.15) to (3.17) naturally associated to the boundary 
hypersurface B. 

Consider the left-side of Eq. (3.75). The boundary condition (3.77) yields ^u^'^las = 
and hence, from Eq. (3.19), 5e"'' = e'^^cr^^^cr'^^las = 0. Thus, we have 

^^'eSB^,Ma^ = Si^'^'eB^.MW (3-78) 

Next consider the right-side of Eq. (3.75). We integrate by parts with respect to the variation 
in the first term to get 

i\^.fjv^:''m = mia^.fcVa^'m - cxs^f/^moise: + c^^y- (3.79) 

Then, using orthogonality relations (3.5) and (3.24), we find that the second term in 
Eq. (3.79) vanishes as follows. First, 

Ct[a^.]Se: = IcisJh - tJsJ = -s^ChhSt' = (3.80) 

since 5t''\Qj: = by the boundary condition (3.77). In addition, 

et[aS.fcK = letMSh; + s^Ss^) = -^Ns^s^Ss^ (3.81) 
Hence, the second term in Eq. (3.79) reduces to 

Ins^s^T^'-'{9)6s'' = (3.82) 

since rj^"'(^) =0. 

Consequently, returning to equation (3.75), we obtain 

SiCB^ie) - m^^s^^eXK^'^mU = -^X^^'i^-^ ^> mei: (3.83) 
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which obviously is satisfied by 

CB^ie) = ^.it^^s^feX^nQ) (3-84) 

and a{^;9,59) — 0. This verifies Proposition 3.1 using the intrinsic Dirichlet boundary 
condition (3.77). 

Finally, from expressions (3.84) for ^'^BJ^O) and (3.63) for Q^^{C,0), we obtain a Hamil- 
tonian (3.76) with the boundary term given by 

HsiCd)^ l^Qai>iC,d)-eB,Jo,se) 

= ^ L ^""^'^V-.r -(^^) - 2t^^s^fXTr{9)). (3.85) 

Hence, the Hamiltonian boundary term takes the form 

HB{^,9) = 8j^^ePciO)dS (3.86) 

where P^{e) = t^sjf'ie) - t^sJlVj'^e) + s^tJlVj'^e). This expression is simphfied by 
the identities (3.3) and (3.14), which yield 

PciO) - f's^a^'^WJ): - Ksy^V.e'', + s^ty^W.el (3.87) 

Thus, we have the following main result. 

Theorem 3.2. For the intrinsic Dirichlet boundary condition (3.77), a Hamiltonian 
conjugate to on E is given by 

H^{C,9)^8jj\{R;{9)-y^Rie))dE + 8j^^ePci&)dS. (3.88) 

On vacuum solutions 9^, the Hamiltonian reduces to the surface integral (3.86), (3.87). 

Note, this Hamiltonian is unique up to adding an arbitrary covector function of the 
Dirichlet boundary data h"' to PJ9). 
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D. Dirichlet Hamiltonian 



On vacuum solutions of the Einstein equations, the Hamiltonian (3.88) with the Dirichlet 
boundary condition (3.77) holding on the timelike hypersurface B bounding a local spacetime 
region V takes a simple form if the frame 9'^ is adapted to the boundary 2-surfaces dT^t and 
e. Let 

< = ta^ K = K^l] = ^a,^ K = ^a^l (3-89) 

which defines a preferred orthonormal frame 'j?^. It follows that 

f = = -51 = s"^?^ = 51 (3.90) 

2 3 

This choice of frame is unique up to rotations of 

Theorem 3.3. For the Dirichlet boundary condition (3.77), the Hamiltonian (3.88) 
conjugate to on evaluated in the orthonormal frame (3.89) for vacuum solutions 
is given by the surface integral 

H''{^;^)^8j^^eP^{^)dS (3.91) 

where 

P^(^) = iV/^V,s„ - ty'^V.s, + s^a''<^V,t,. (3.92) 

We refer to H^{^', 'd) as the Dirichlet Hamiltonian for the gravitational field in the local 
spacetime region V, and P^i'd) as the Dirichlet symplectic vector associated to the boundary 
2-surfaces dT,t. Note that, since ^" lies in B, only the first two terms in P^{'&) contribute 
to H^{^;-d). The significance of the full expression for P^('f?) and its resulting geometrical 
properties are discussed in Ref. [8]. 

The general form of the Hamiltonian boundary term (3.86) and (3.87) differs from the 
special form (3.91) and (3.92) when evaluated in an orthonormal frame not adapted to dJ^t 
and In particular, we obtain the relation 
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PciO) - P^i^) - ^'t^d,s^ + t/^d,s^ - s/^d,t^ (3.93) 

and so the general form of the symplectic vector P Jy&) differs from P^{'&) by various gradient 
terms. These terms can be understood by considering a change of orthonormal frame 

- U'X (3.94) 

where U'^^ is a SO (3,1) transformation acting in the frame bundle of the spacetime {M,g) 
at dTit- Such transformations are defined by U = U^pTj^^r]'^^ and det{U) = 1, where 
C/~^^„ is the inverse of U'^„, given by U"^'^,,!/'^^ — St — U'^,.U~^'^^. 

The transformations (3.94) arc a gauge symmetry of the tetrad formulation for General 
Relativity. Under the change of orthonormal frame, one has 

(3.95) 

and so, through substitution of the transformation (3.95) into the curvature (3.51), 

Ra,:{0) - Ra,:m = U-'^U'^.Rjie) (3.96) 

after cancellations of terms. Hence the Lagrangian (3.52) for the field variable is gauge 
invariant. As a consequence, it is straightforward to see that the symplectic structure given 
by the symplectic potential (3.55) and current (3.66) must be gauge invariant. In particular, 
note that one has SV^^^Ue) = U'^'^^U'' p5V^J{e) where the gradient term from Eq. (3.95) 
drops out of the variation since it has no dependence on This explicitly establishes the 
gauge invariance of ^ahc^^^SO) and hence of uj^^^{6,5i6,5-26). 

However, the Noether charge (3.63) fails to be gauge invariant due to its explicit de- 
pendence on the frame connection. Consequently, it follows that the gradient terms in the 
symplectic vector (3.93) originate from a gauge transformation on the frame connection un- 
der (3.94) as given by a transformation relating the adapted orthonormal frame to a general 
orthonormal frame, 9^, = U^r,t„ + U^.a^^ + C/^cr^^ + U^.s. 

'a U a lo 2 a 4a 
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The gauge invariance of the symplectic structure arising from the tetrad formulation of 
the Lagrangian means that the symplectic potential (3.55) and current (3.66) are equivalent 
to manifestly gauge-invariant expressions derived using the metric formulation of General 
Relativity with gi^^ as the field variable. It can be shown that one has 

e>a6c(^, ^9) = e>a6c(^, + 69[,V'.,] , 5g) (3.97) 

where V'ablfl') ^o) is a locally constructed 2-form, and so the symplectic potentials are equiva- 
lent to within an exact 3-form. This contributes a boundary term to the presymplectic form 
obtained from the metric Lagrangian L^j^^^[g) = ^abcdi9)-^i9)^ 

n,:{g,dg,C^g)=nj:{9,d9,CS+ [ e''\d^Jg,C^g)-C^i;Jg,dg))dS. (3.98) 

Correspondingly, the Noether charge 2-form Q^bi^^ q) arising in the metric formulation differs 
from QabiC^ ^) ill the tetrad formulation by the term 2ip^f^{g, C^g) . Explicitly, using the metric 
Lagrangian, one finds that [11] 

\e'Qa,{i. 9)- -4t[c5.]^Vr = 4eV^V,5, - 2{s''C^t^ + f (3.99) 

Here the first term in Eq. (3.99) is simply the Noether charge (3.63) evaluated in the adapted 
orthonormal frame (3.89), 

^eTr{^y'e^^^{g)^y: = i6i\s^^"''v^r, = leet'^v^s, (3.100) 

since V^i^'V^^s^-^) = by hypersurface orthogonality of s^. The second term in Eq. (3.99) 
simplifies through the hypersurface orthogonality relations (3.6) and (3.7), leading to 

fC^s^ = m'^V^^s^^ + a,(e%)) = (3.101) 

and 

N 

s'^C^t^ = -e-dal^ (3.102) 

where a, (3,N are scalar functions defined by 
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s = ads, t = -N(dt + (3ds) 



(3.103) 



with 



C^t = ed^t = 1, = f = 0. 



(3.104) 



Hence we obtain the relation 



(3.105) 



A similar relation can be shown to hold between the respective symplectic vectors arising 
in the tetrad and metric Hamiltonian formulations of General Relativity. In particular, by 
direct calculation with g^^ as the field variable, one finds that the full Dirichlet boundary 
condition (3.34) yields a Hamiltonian conjugate to ^" on S whose boundary term is given 



This differs from the symplectic vector in the tetrad formulation by the same gradient 
term occurring in the Noether charges (3.105). The extrinsic part (3.38) of the Dirichlet 
boundary condition is necessary in obtaining this Hamiltonian, because of the boundary 
term in the presymplectic form (3.98). Interestingly, in the case when is orthogonal to 
Ej, then /? = 0, and one finds that the weaker, intrinsic Dirichlet boundary condition (3.37) 
is sufficient for existence of the metric Hamiltonian (3.106) and (3.107). Moreover, in this 
case the presymplectic form (3.98) and symplectic vector (3.107) are exactly the same as 
those obtained in the tetrad formulation using the adapted orthonormal frame (3.89). 

An expression for the Dirichlet Hamiltonian boundary term (3.106) in terms of the 
standard ADM canonical variables associated to S, and its relation to quasilocal quantities 
of Brown and York [16,17], will be derived in Sec. IIIF. 



by 




(3.106) 



where 



(3.107) 
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E. Determination of allowed boundary conditions 

A necessary and sufficient condition [10] on variations 69'^ for existence of a Hamiltonian 
(3.76) conjugate to on E is given by the antisymmetrized variation of the equation (3.72) 
for the boundary term. This yields 

e'^'Cu^abci^, Sie, 620)19^ = ^A^"(e, 0; 6,6, 52^)|aE (3.108) 

with 

/3"(e, 9; 6,9, 629) = e^Xa.i^, 9; 6^9) - e'''62a,{i, 9; 6,9). (3.109) 

To begin, we simphfy the expression (3.66) for u!^f^^{9,6i9,629). First, using Eq. (3.74) for 
®a6c(^' taking into account the orthogonahty ^"s^ = 0, we have 

^^""'eeaJO, 69)^ ehK6{s^h'^^VJl) (3.110) 

through the frame decomposition (3.41) and the relation r/''(^) = 9"^^^} = -9"''^^} ^9^^. 
Now we substitute the identity Jf^ — h^'^^\/^9l — s^s'^V J)^^^ and then use the relations 
h"'^6s^ = 0, fi'^s^ = to simplify the term 

h;Si\\h'^^^'^cdl) = 0- (3.111) 

Thus, Eq. (3.110) becomes 

^e"'r0„.c(^, = c\h'"'s{s^h;h:^v^9i). (3.112) 

Finally, we substitute ^'^t^ = ^^'^''C^abd^)- 
Hence, we obtain 

rsQabci^, 69) = 8e^,^{h)h;6K/ (3.113) 

where we define 

Ka' = «>'"OVe^; = s^h:rr{9). (3.114) 
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Note that s"'Kj^ — and s^Kj^ — 0. Prom Eq. (3.113), by taking an antisymmetric 
variation and then using Eq. (3.46) for the variation of e^(,g(/i), we have 

= Kdh) {{Sih; - h;h:5,h;)5,K/ - {5,h; - hX^^Ky.Kj-) . (3.115) 

Substitution of this expression into equation (3.108) yields the following result. 

Theorem 3.4. A Hamiltonian conjugate to on E exists for variations 50^ with support 
on if and only if 

for some ^i^) ^26') of the form (3.109). The Hamiltonian is given by 

H^{i;0) = 8 - yMm^ + HB{i;9) (3.117) 
with boundary term 

Hb{U) = lQa,iU)-eB,,M (3-118) 
where B^^^{6) is determined from the equation 

CiVsSB^je) - Se^^^mX/) = Vsd^,a,^{^, 6- 66). (3.119) 

Thus, equation (3.116) determines the allowed boundary conditions on variations 69'^ 
for existence of a Hamiltonian formulation (3.117) for the vacuum Einstein equations. To 

proceed, we now parallel the analysis of the similar boundary condition determining equation 
for the Maxwell equations in Sec. II C. 

Two obvious solutions of the determining equation (3.116) with ^""{^,9; 5i9,620) — 
are given by 5/i^"|aEt = 0,t > 0, which is the Dirichlet boundary condition (3.77) already 
considered; and by 
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5i^/|as, = 0, t>0 (3.120) 

which we call the Neumann boundary condition. For the boundary condition (3.120), it 
follows from Eqs. (3.76) and (3.119) that the corresponding Hamiltonian boundary term is 
given by 

H''{C;e)^8j^J'P^{e)dS (3.121) 

where 

Pci&) - t^^u^nS) - tX^'^X (3-122) 

by a derivation similar to Eq. (3.87). In the orthonormal frame (3.89) adapted to the 
boundary 2-surfaces dTif, we have 

P,^(^) = t"^V,s„ (3.123) 

We refer to this as the Neumann symplectic vector associated to the boundary 2-surfaces 
dT,t- Moreover, in this frame the Neumann boundary condition (3.120) becomes 

SKf\ej:, = 5{h'''KjU, = 0, t>0 (3.124) 

in terms of 

Ka, = KX'^cSa (3.125) 

which is the extrinsic curvature of the timelike boundary hypersurface B in {M,g). Thus, 
geometrically, the Neumann boundary condition corresponds to fixing the frame components 
of the boundary hypersurface extrinsic curvature, 

5(/i'"'/i„"^V,sJ = on B. (3.126) 

These components measure the rotation and boost of the hypersurface normal with respect 
to the frame hj^ under displacement on B. 
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We now investigate more general boundary conditions. Note that, on the left-side of the 
determining equation (3.116), e"''^^'^cu^f^^{9,Si9,S29) involves only the field variations Vs^h^ 
and VbSKJ^ = VbS^s^^H^T J^^ {9)). We call and Kj^ the symplectic boundary data at 
dTif and consider boundary conditions of the form 

5T:{h:,K:)u^o,t>o (s.m) 

where ^^i^c^ ^c) locally constructed as an algebraic expression in terms of the symplec- 
tic boundary data and fixed quantities (including the spacetime coordinates). We call (3.127) 
a mixed Dirichlet-Neumann boundary condition if ^a(^/, K^) is a constant-coefficient lin- 
ear combination of the parts VashJ^, VthJ^, Vq-eKJ^, VtKj^ of the Dirichlet and Neumann 
boundary data in (3.77) and (3.120). Here the projections with respect to Vqj: and Vt remove 
all components proportional to s^. 

An analysis of the boundary condition determining equation (3.116), given later, leads 
to the following main results. 

Theorem 3.5. The only allowed mixed Dirichlet-Neumann boundary conditions for 
existence of a Hamiltonian (3.76) conjugate to on E are given by 

VBiaoSKj" + boShJ')U, = 0, t>0, (3.128) 

or equivalently 

niK^K:)^a,K:+b,h: (3.129) 

for constants Oq, &o (and P"' — in Eq. (3.116)). In the cases Qq — or bg — 0, respectively 
Dirichlet or Neumann boundary conditions, the corresponding Hamiltonian boundary terms 
are given by Eqs. (3.91) and (3.87), and Eqs. (3.121) and (3.122). In the case ao 7^ 0, 6o 7^ 0, 
the corresponding Hamiltonian boundary term is given by 

H^^{C;9)^8f CPcWS (3.130) 

where, now, 

39 



^cW = ^fW + 6V. (3.131) 
do 

(Note, the boundary terms here are unique up to adding an arbitrary covector function of 
the boundary data (3.129) to P^{9).) 

A similar covariant derivation of the pure Dirichlet and Neumann boundary terms is 
presented in Ref. [18,19] from a different perspective. In Theorem 3.5, note that Eq. (3.128) 
represents a one-parameter ao/&o family of boundary conditions. In particular, in contrast 
to the two-parameter family of analogous mixed Dirichlet- Neumann boundary conditions 
(2.52) allowed for the Maxwell equations, here decompositions of the symplectic boundary 
data with respect to Vdn and do not yield boundary conditions satisfying the determining 
equation (3.116). 

The form of the mixed Dirichlet-Neumann boundary condition (3.128) suggests we also 
consider boundary conditions specified by a trace part and trace-free part with respect to 
the boundary hypersurface frame hj^: 

S:F{h:,K:)\aj:,^0, <5^^(/.;, ir/)|aE. = 0, t>0 (3.132) 

with h^P^{h^ ,K^) ~ 0. Taking the trace of the symplectic boundary data variations 
yields 

K^^a =^K + Wii'/^V (3-133) 

and 

h;5h: = 5hi\h\ (3.134) 
where K = h "^K ^ is the trace of and h = det(h ^) is the determinant of the frame 

a a V a / 

components hj^ in local coordinates. 

Theorem 3.6. Allowed boundary conditions (3.132) for the existence of a Hamiltonian 
(3.76) conjugate to ^" on E are given by 
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VB5{Kj'-^-h^^K)\a^, = Q, t>0, (3.135) 

{aoSK + boS\n\h\)\9j:,^0, t > 0, (3.136) 

or equivalently 

T{h, K) = aoK + 6o In P^^K , K;) = VsiKj- - ^h^K) (3.137) 

for constants ao,&o (and /3" = m (3.116)). The corresponding Hamiltonian boundary 
term is given by 

Hb{^: 9)^s[ edPciO) + lPci(^) + "^-^MS (3.138) 

O O 0,Q 

( which is unique up to adding a term depending on an arbitrary covector function of the 
boundary data (3.137)). 

Finally, we remark that the mixed boundary conditions in Theorems 3.5 and 3.6 admit 
the following two generalizations. 
First, 

^'aiK^ = a{x, K, In IkDrsiKj" - ^-hJ^K), (3.139) 
P(h,K) ^b{x,K, In \h\), d^by^O (3.140) 

for arbitrary functions a{x, K,ln\h\), b{x, K,ln\h\). 
Second, 

KiK^ = a{x, K, In \h\)VB{K^'' + b{x, K, In |/i|)/i/), b + -\k (3.141) 

for an arbitrary function 6(x, In |/i|), with the function a(a;, In |/i|) now satisfying the 
linear partial differential equation 

obtained from the determining equation (3.116). The general form of a is given by solving 
the characteristic ordinary differential equations 
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d\n\h\ dK da , , 
^ = = — 3.143 

V ^ + ^in|.|^ + 3d ^b)a 

in terms of the variables In \h\, K, a. For instance, if b is taken to be linear homogeneous in 

2 A 

K , then one has b = XK , a = f{x, XK In ^ ^^-i ^ where A = const and / is an arbitrary 

function. 

Proofs of Theorems: 

Since any boundary condition locally constructed from the symplectic boundary data is 
hnear homogeneous in Vs^hJ^ and VbSKJ^, we begin by finding all such solutions of the 
determining equation (3.116). 

First we show that = 0. The right-side of Eq. (3.116) necessarily involves terms with 
at least one derivative on 59^, while only first-order derivatives of 69'^ appear on the left-side 
of Eq. (3.116) through 

VbSK/ = s^h^Srrie) (3.144) 
due to Eqs. (3.114) and (3.47). Thus, for a balance in numbers of derivatives, we must have 

«a = O^aU^M (3-145) 

for some Oi^^{9) locally constructed out of 9'^ and fixed quantities. This yields, for the 
antisymmetrized variation of cc^, 

(^a = <^a',.mQ'Mc (3-146) 

where Q;^^^j^(^) = —c^a^ni^) is the curl of oc^JJ)) with respect to 9"^. Then, using Eqs. (3.109) 
and (3.50), we collect all the terms on the left-side of Eq. (3.116) linearly independent of 
ShJ^ , SKJ^. Through Eqs. (3.144) and (3.49), the coefficients of these terms yield 



«c«7^'^(^) = 0, «,<l^(^)-0, s^a^'X{9)^0 (3.147) 

where 
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= = (3.148) 

These algebraic equations are straightforward to solve, leading to 

= t^'c^itm + a^tM + t'fa.lM + t^'o.4M (3.149) 

for some 

ca c a bd abc a b c edf a a e ca c a bd 

(3.150) 

Then, returning to Eq. (3.145), we note that Q;^^(^) = ^a^^ct^Q for some that is locally 
constructed only from fixed quantities since it is a scalar expression. Thus we immediately 
obtain 

abc ^ eye'^'^ia^ -a^ - / + S^). (3.151) 
By equating expressions (3.149) and (3.151), we find that after some algebraic analysis 

< = €«o, (3.152) 

for some a^^. Now, substitution of expression (3.152) back into Eqs. (3.149) and (3.151) 
easily leads to the result 

ca abc a ca r\ /o-iroA 

= ^^t^v = ^^i^v = «4^^ = 0. (3.153) 

Hence, from Eq. (3.149), we have oi^^y{0) = and so (3^ = 0. This establishes that P"' = 0. 
Consequently, the determining equation (3.116) reduces to 

{5,h; - h;h:5,h;)5,K/ - {5,h; - h;h:5,h:)5,Kj^ = o, (3.154) 

which is equivalent to 

K\{i^iK^2K,^ - S^h^S.K/) = 0. (3.155) 
Then the algebraic solution of Eq. (3.155) in terms of Ve^h^ and Vb^KJ^ has the form 
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^iZPb^K + ^2ZVb5K,^ = (3.156) 
for some coefficient tensors 11'^'' such that 

Ktha%{ = Ka\:K\- (3-157) 

It is straightforward to show that Eq. (3.157) holds iff 

K = K-KK + K'K (3-158) 

where fi^j^h^h^ — ^'^hj^hj' is the symmetric part of U.'^l in the index pairs (a, //) and 
(6, v), and where fl|^ = /i^^H^^, = hj^^t^al is the trace part of 11^^ in the frame hj^. Thus, 
we have established the following result. 

Lemma 3.7. All solutions of the determining equation (3.116) for allowed boundary 
conditions that are linear homogeneous in Vs^hJ^ and Vb^KJ^ have the form (3.156) where 
the coefficient tensors are given by Eq. (3.158). 

Now, for mixed Dirichlet-Neumann boundary conditions, we take 

ni^^ = ai(PaE)„V/ + ao{V^)Xt' (3-159) 

Ti^t = h{Vo^)X: + ho{Vt)Xe (3.160) 
where ao, Oi, 6o, &i are constants. Then the requirement (3.158) leads directly to 

ao = Oi, &o — ^1- (3.161) 

Substitution of Eqs. (3.159) to (3.161) into Eq. (3.156) yields the mixed Dirichlet-Neumann 
boundary conditions (3.128). 

Finally, the other boundary conditions (3.135) and (3.136) arise from 

^.-I^K'K'-IkKK'. (3-162) 

= hX', (3-163) 

which are easily verified to satisfy the requirement (3.158). 
This completes the proofs of Theorems 3.5 and 3.6. □ 

As a concluding remark, wc note that Lemma 3.7 yields the following necessary and 
sufficient determining equations for finding all boundary conditions (3.127). 
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Lemma 3.8. All allowed boundary conditions of the form ^^(/i/, i^/)|as = for 

existence of a Hamiltonian conjugate to on S are given by the solutions of the equations 



— ^ u^h,,^ = h,;hj, (3.164) 

= K'hJ- (3-165) 

We will leave a general analysis of the boundary condition determining equations (3.164) 
and (3.165) for elsewhere. 



F. Relation between covariant and canonical Hamiltonians and boundary terms 

To conclude this section, we first give a brief discussion of the Hamiltonian field equations 
for General Relativity using the covariant symplectic structure and Noether charge Hamil- 
tonian in Sec. HIC. Then we discuss the boundary terms in the Dirichlet and Neumann 
Hamiltonians expressed in standard ADM canonical variables. 

For a Hamiltonian (3.76) conjugate to on E, the associated field equations are obtained 
through the presymplectic form (3.69) by the variational principle 

Qe(^, 59, c^e) - H'^ii- e, 56) = 1^ 8^%(R;(e) - ^e;)5e'^ = o (3.166) 

for arbitrary variations S9'^\j:. These field equations split into evolution equations and con- 
straint equations with respect to E corresponding to a decomposition of 9'^ into dynamical 
and non-dynamical components determined by [4] the degeneracy of 

nJ:{9,^^9,^29) = l^u;^J9, 6^9,629) - ^e'''"'n^u;^,^{9,dAd29)dE. (3.167) 

For this purpose, it is convenient to partially fix the SO (3,1) local gauge freedom in 9'^ 
analogously to conditions (3.15) to (3.17) by choosing the frame components n'^ — n"'9'^ 
to be fixed constants on M. Then, through a simplification of the symplectic current here 
similar to Eq. (3.115), we obtain 
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^e"^"^n^u;„^^(^, S^O, 626) = - q^^q^^ 5iqJ)52Pd + (^2?/ - g/C'^2g/)5ip/ (3.168) 

where we define 

with q^^ = g^^ + nji^, which are counterparts of hj^ , Kj^ associated to the spacehke hy- 
persurface E. (Geometrically, qj^ is a frame for the hypersurface metric q^^ — qj^q^'n^vi 
while pj^ represents the frame components of the hypersurface extrinsic curvature K,^^ — 
^aPb'^ixv ~ ^a^^ J^b ) Hcncc, the degeneracy directions of the presymplectic form (3.167) 
are given by variations 56^ satisfying 

%a7kl) = ^Pa" = 0. (3.170) 

where |g| = det(qj^) and 51n|g| = q^bq^ . This immediately leads to the following result, 
analogous to the discussion in Sec. II A for the Maxwell equations. 

Proposition 3.9. For a Hamiltonian H-s{^;9) conjugate to on E, let ^kero; denote 
the field components of 9^ and Tj^^{9) invariant under the symplectic degeneracy directions 
(3.170), and let denote the remaining field components modulo 6]^eTuj- Then there is 
corresponding decomposition of the Hamiltonian field equations given by H^{^; 9, 59keru>) — 
and H^{^■,9,59^J) = flj:{9,59^, C^9), which respectively yield 

nRj'{9) - ^n^R{9) = 0, (3.171) 
qaKi^) = 0. (3.172) 

These field equations arise equivalently by variation of Vn9^ = —n^n^ and Vt,9^ = qj^ in 
the Lagrangian (3.52). 

Equations (3.171) and (3.172) are the frame components of the standard 3+1 split of 
the vacuum Einstein equations [12] into constraint equations and time-evolution equations 
for the hypersurface metric q^f^. Thus, the Hamiltonian field equations given by the varia- 
tional principle (3.166) constitute a covariant formulation of the standard ADM Hamiltonian 
equations for General Relativity. 
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With respect to the spacehke hypersurface E, one has a decomposition of ^" into normal 

and tangential parts 

C^Nn+N" (3.173) 

where jV" = 'Ps(C") ^'^'^ — ~C'^a define the lapse and shift of the time fiow vector field 
By use of the Gauss-Codacci equations, we straightforwardly see that the volume part 
of a Hamiltonian (3.76) conjugate to on E is given by the "pure constraint form" [11] 

H{i- e) = A j M{n + IC^- /C„^/C"') + 2Af'{v'lC^^ - VJC) dL (3.174) 

where TZ^^ and K,^^ are the Ricci curvature and extrinsic curvature of the metric Vj^Qab — lab^ 
TZ = TZ^ and K, = K,^ are the corresponding scalar curvatures, and is the derivative 
operator associated with q^. (An analogous result holds more generally for any diffeo- 
morphism covariant Lagrangian field theory [4]. ) This demonstrates, explicitly, that our 
covariant analysis of allowed boundary conditions and corresponding boundary terms for 
General Relativity in Sees. Ill B and III E is equivalent to a canonical analysis of the ADM 
Hamiltonian. 

Now, consider Dirichlet or Neumann boundary conditions imposed at the 2-surfaces SE^, 
for i > 0. On solutions of the Hamiltonian field equations, the total Hamiltonian H-s{^;9) 
reduces, respectively, to the Dirichlet and Neumann boundary terms (3.91) and (3.121). Let 
u"" denote the outward unit normal to dT,t in ^t- Let i?!^ be an orthonormal frame adapted 
to Et given by 

= ^a, K = ^a, = ^a,^ K = ^a^l (3-175) 

which is related to the frame adapted to i3 by a boost in the normal space T-'-(9Et) to 
the boundary 2-surface 9Et, 

= cosh X + ti" sinh X, = ti" cosh x + sinh (3.176) 

Through the corresponding boost relation (3.95) applied to the symplectic vectors (3.92) 
and (3.123), the Hamiltonian boundary terms take the respective form 
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eP^i^) = e{n<7f^V,u^ - n^'^V^u, + u^'^V.n, - af^V,x), (3-177) 
rP^(^) = r(n"^V,«„-W,x)- (3.178) 

These expressions can be simplified in terms of the hypersurface metric q^^^, extrinsic cm^va- 
ture /C^^, and acceleration Oj, = n^^Vjif^. We find that 

CP^m =Nk- Af\%, - Xp^x, (3.179) 
fPf (^) = -J^u% - N'-u'lC^, - Nd,x - M^d^x - Mld^X. (3.180) 

where jVy = Pas (A/""), N — uj^f"' are the tangential and normal parts of the shift with 
respect to dllf> and k — c^^^V^Wf, is the mean extrinsic curvature of dUt in '^t- 

We note that this form of the Dirichlet and Neumann boundary terms (3.179) and (3.180) 
agrees with the canonical analysis of boundary terms for the ADM Hamiltonian carried out 
in Refs. [16,17]. Moreover, in the case when is orthogonal to B, i.e. x = 0, the surface 
integral Jqy^ ^'^P^{^)dS for suitable choice of ^'^ reproduces Brown and York's expressions for 
quasilocal energy, normal momentum, and tangential momentum quantities (respectively, 

= = 0; = 0,A^jJ = 0,Ar^ = 1; = = 0,A^jJ 7^ 0). Further discussion 

of quasilocal quantities associated to the Dirichlet and Neumann symplectic vectors (3.92) 
and (3.123) will be left for elsewhere. 

IV. CONCLUDING REMARKS 

In this paper we have given a mathematical investigation of boundary conditions on the 
gravitational field required for the existence of a well-defined covariant Hamiltonian varia- 
tional principle for General Relativity when spatial boundaries are considered, with a fixed 
time-flow vector field. In particular, a main result is that we obtain a covariant derivation of 
Dirichlet, Neumann, and mixed type boundary conditions for the gravitational field in any 
fixed spatially bounded region of spacetime. We show that the resulting Dirichlet and Neu- 
mann Hamiltonians lead to covariant Hamiltonian field equations which are equivalent to the 
standard 3-1-1 split of the Einstein equations into constraint equations and time-evolution 
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equations. In addition, we obtain a uniqueness result for the allowed boundary conditions 
based on the covariant symplectic structure associated to the Einstein equations. 

However, we do not address the purely analytical issue of whether the boundary-initial 
value problem for the Einstein equations is well-posed with these boundary conditions (i.e. 
do there exist solutions of the Einstein equations satisfying the boundary conditions, initial 
conditions, and constraints). For work in that direction, see e.g. Ref. [20]. 

A further interesting generalization of our results would be to treat a spacetime region 
whose spatial boundary is dynamical e.g. a black- hole horizon or Cauchy boundary. We 
note that boundary conditions for this situation may be investigated by allowing the time- 
flow vector field to depend on the spacetime metric instead of being a fixed quantity. This 
analysis will be pursued elsewhere. 
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First of all, consider in n spacetime dimensions a general Lagrangian field theory for a 
set of fields denoted collectively by (/>. It will be assumed that these fields are defined as 
sections of a vector bundle E over the spacetime manifold M, using local coordinates on M 
and E. The theory will be assumed to be described by a Lagrangian n-form L(0) that is 
locally constructed out of the fields and their partial derivatives d^cf) to some finite order 
k (and fixed background structure, if any, on M and E) . 

The Lagrangian L{(f)) provides a variational principle 



which yields the field equations = obtained as the stationary points of S'(0), 



APPENDIX: NOETHER CHARGE METHOD 




(Al) 
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dS{(l)) = / = / £{(l))S(f) = (A2) 

Jm Jm 

under variations 50 of with compact support on M. For arbitrary variations 50, which 
are not restricted to have compact support, one then has a variational identity 

L'{(j), 6(1)) = 6L{(j)) = £{(t))5(t) + (ie(0, 50) (A3) 

where 0(0,50) is an {n — l)-form, caUed the symplectic potential, derived through formal 
integration by parts. This yields a well-defined locally constructed formula for f (0) and 
0(0, 50) in terms of 0, 50, and their partial derivatives to a finite order. The symplectic 
potential is used to define the presymplectic form on a fixed hypersurface E, 

f^E(0, 5i0, 520) = ^ a;(0, 5i0, 520) (A4) 

in terms of the symplectic current {n — l)-form uj given by 

u;(0, 5i0, 520) = 5ie(0, 520) - 520(0, 5i0). (A5) 

The symplectic current satisfies ('icj(0, 5i0, 520) = £^'(0, 520)5i0 — £'(0, 5i0)520 with 
E'{(f), 50) = SS{4>). Hence, a;(0, 5i0, 520) is closed for variations on solutions, 

ciu;($,5i*,52$) = (A6) 

where $ denotes restricted to satisfy ^(0) = 0, and 5$ denotes 50 restricted to satisfy 
£^'($, 50) = 0, i.e. 5$ is, formally, a tangent vector field on the space of solutions. Con- 
sequently, f^E(0, 5i0, 520) is unchanged by deformations of the spacelike surface S in any 
compact region of M. 

In the previous constructions, a change of coordinates on M or £■ leaves £^(0) unchanged 
[22], while 0(0, 50) changes in general by an exact locally constructed (n — l)-form (ii/(0, 50). 
However, one can show that [4] if the Lagrangian 1/(0) is at most second order in partial 
derivatives d'^cf) {k < 2) of 0, then 0(0, 50) is independent of the choice of coordinates 
on M and E and thus the presymplectic form ns(0, 5i0, 520) is then coordinate invariant. 
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Moreover, note that L(0) can be freely changed by addition of a locally constructed exact 
form dfi{(j)), without affecting the field equations S{(f)). This changes 6(0, (50) by addition of 
a locally constructed (n — l)-form ddiJ,{4>), but leaves 0^(0, Si4>, ^20) unchanged. Therefore, 
up to its dependence on E, the symplectic structure Jls(0, 5i0, 52(p) is uniquely determined 
by £{4>) in this situation. 

Now consider a complete, nowhere vanishing vector field ^ on M. It will be assumed 
that there exists a well-defined Lie derivative acting on associated to the diffeomorphism 
generated by $, on M. Let E be a connected region contained in a fixed hypersurface in M 
with a closed boundary 9E. (Note, if E is simply-connected, 9E is a closed n — 2-surface 
in M bounding E. If E is multiply-connected, then 9E is a disjoint union of closed n — 2- 
surfaces. Also, if E extends to "infinity", then 9E contains a corresponding "asymptotic 
boundary" n — 2-surface.) 

Definition A.l. A Hamiltonian conjugate to ^ on E is a function ifs(^;0) = /s^(^;0) 
for some locally constructed (n — l)-form Tt{i; 4>) such that, on solutions 



where denotes the Lie derivative, and H^{^; 4>, 64>) = 6H-s{^; 4>). 

This is a covariant formulation of the standard Hamiltonian symplectic structure, with 
the "time" direction defined by ^, called the time fiow vector field. (In particular, Ref. 
[4] outlines a construction of a standard phase space and Hamiltonian equations of motion 
determined from this covariant structure.) Note that a Hamiltonian, if it exists, is auto- 
matically conserved along ^ for solutions $, i.e. £^Hj^{^; $) = 0. It will now be shown that 
when ^ is a "time symmetry" of the Lagrangian L{4>), then a Hamiltonian conjugate to ^ 
exists and is simply given by the Noether charge associated to ^. 

Given any vector field ( on M, consider the variation 6^(j) = ^c.'P- If this is a symmetry 
of the Lagrangian, so that 



i/^(e;$,50) = ^^E($,50,A^) 



(AT) 



(A8) 
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by means of the identity >C^L(0) = d{i(^L{(l))), then one can define a conserved Noether 
current {n — l)-form J{(; (p) by 

J(C;0) = e(0,£c0)-^c^(0) (A9) 

where is the interior product. Conservation of this current simply means that, on solutions 
J(C; <p) is closed 

dJ((; $) = de($, - d{i<;L{^)) = L'($, 5$) - CcH^) = (AlO) 

through the symmetry condition (A8). The integral of J(C; ^) over E defines the Noether 
charge 

QdO- (All) 

One finds that the "time" derivative of this charge with respect to ^ is given by 

AQe(C)= / A-^(C;^)= / kdJ{C;^) + d{^J{C;^))^ { hJ{C;^) (A12) 

where i^J{C', <P) is called the /?ua; of the Noether current. Hence, if the flux vanishes on 9E, 
then the charge is conserved for solutions 

Examples of field theories which admit a symmetry 5^0 = £^0 are (i) any generally- 
covariant theory on a fixed, background spacetime (M, g) with an isometry vector field ( 
(i.e. jC(^g — 0), where is purely a function of g, (f) and its metric-covariant derivatives 
V0; (ii) any diffeomorphism-covariant theory, whose field variables (p include the spacetime 
metric g, where L{(f)) is purely a function of 0, curvature tensor of g, and their metric- 
covariant derivatives. 

For a diffeomorphism-covariant theory, 5<^0 is a symmetry for all vector fields Conse- 
quently, since J{(; 0) is locally constructed out of (, one can show that in this case [4] 

J(C;*) = c^g(C;$) (A13) 

for some locally constructed (n — 2)-form Q{C; 0), called the Noether current potential. Then 
the Noether charge reduces to a surface integral 
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QdO- [ J{C-A)=i Q(C;0)- (Ai4) 

In contrast, for a generally-covariant theory, J(C; is related to the conserved stress- 
energy tensor T(0) defined by considering variations of 

*6gL{<P) = -^T{<P)5g + *de{<P, dg). (A15) 

One can show that [9], on solutions $, 

J(C;$) = *^c^($) + rfr(C;$) (A16) 

for some locally constructed (n — 2)-form r(C; 0). 

Proposition A. 2. For any symmetry 5^(1) — C(^(f) admitted by a Lagrangian L{(f)), the 
field equations and symplectic potential satisfy 

6<:S{<i)) = C^m, (A17) 

5ce(0, 50) = £ce>(0, + #(C; 0, (ai8) 

where 0, (^0) «5 some locally constructed {n — 2)-form. 
Proof : 

Consider an arbitrary variation of the Lagrangian symmetry condition ( A8) , 

= 5{5(^L{(i)) - C^L{(t))) = 5^5L{(t)) - C^SL{(I)). (A19) 
From Eq. (A3), one has 

£^(5L(0)) = £^(£(0)50) + C(;dQ{(j), 6(j)) = {C^S{(j)))6(j) + S{(j))6C(;(i) + ^£^9(0, 50), (A20) 
and similarly 

d^{5L{(P)) = {5^S{(P))d(P + £{(P)dC^(P + d6^Q{(P, 50), (A21) 
since 5^0 = £^0. Hence, Eq. (A19) yields 
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{SiS{(l)) - Ci£{(f)))S(f) = ci(Ae(0, 5(f)) - S^e{(j), S(j))) (A22) 

holding for all 50. By taking 50 to have compact support and integrating the equation 
(A22) over M, one obtains j^{5c^£{(t)) — Cc^8{(t)))5(p — which immediately yields Eq. (A17). 
Then Eq. (A22) shows that C^Q{4>, — 5(^Q{4>, is a closed {n — l)-form holding for all 
(f). Since this expression is locally constructed in terms of 0, it follows that [21,22] Eq. (A18) 
holds. n 

From these results, one finds that the variation of the Noether current is given by 

J'(C; 0, 64>) = dJ{C, 4>)= 5e(0, £c0) - ^i^m 

= a;(0, 50, C^(t>) + 5^6(0, 50) - ic(de(0, 50) + £:(0)50) 

= ^(0, 50, £c0) - ^d^m4>) + c^(^ce(0, 50) + V^(C; 0, 50)) (A23) 

using the identity i,^{dQ{(p, 6(p)) = £^6(0, 50) — d{i(^Q{(j), 6(j))). 
Lemma A. 3. On solutions 

Qs($, A^, Scj>)^- [ J'(C; 50) + / zce($, 50) + V(C; ^, <^0)- (A24) 

T/iu5, /or variations 5(1) with compact support in the interior ofE, i.e. 50|gs = 0, 

Cc^, 5(t>) = - J'(C; 50). (A25) 

One can then apply this result to the time flow vector field ^ = ^ to obtain a Hamiltonian. 

Theorem A. 4. The Noether current J{C',(f>) yields a Hamiltonian conjugate to ^ on 

S given by /Js('C;0) = under compact support variations 5(p. For solutions 

-f^s(Cj^) — ^he conserved Noether charge associated to ^. 

For variations 50 without compact support, there exists a Hamiltonian if and only if one 
can find a locally constructed {n — 2)-form B{^;(f)) such that 

/ B'iC, 50) - i^e($, 50) - V'lC; S(f>) = (A26) 
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where B'{^; 0, 5(f)) = SB{^;(f)). If one restricts to variations 5(f) — 5$, then by considering a 
second variation and antisymmetrizing in this equation, one obtains the necessary condition 

/ 5i(%e($,52$) + ^/'(e;*,52$))-52(%e($,5i$) + ^/;(e;$,5i^')) = (A27) 

for existence of B{^; (f>). This condition can also be shown to be sufficient [10]. 

Definition A. 5. An allowed boundary condition on (f) is a set of field components J^(0)|as 
locally constructed from (f>, partial derivatives d^(f), and spacetime quantities associated to 
such that for all variations 5(f> satisfying J-''{(f>,5(f>)\dj: = 0, where J-''{(f>,5(f>) = 
5T{(f)), there exists a Hamiltonian Hs{^;(f)) conjugate to ^ on E. 

One now has the following main result. 

Theorem A. 6. A Hamiltonian conjugate to ^ on E exists under variations 5(f) without 
compact support if and only if 

(f i.uj(^,5i^,52^) = <f V^'fC; ^2$) (A28) 

Jdi: Jdll 

on solutions where ifj\^; 0, ^20) = ^iV^lCi 4'i ^20) — ^2V'(Cj 'Pi Si4>)- This determines the 
allowed boundary conditions J-{(f>)\ai: for the field equations to admit a covariant Hamiltonian 
formulation. Then the Hamiltonian is 

Hj:(^; (f)) = 1^ J(e; 0) - dB{^; 0) (A29) 

with B{^;(f)) given by Eq. (A26) up to an arbitrary function of the boundary data J-'{(f)) and 

^. Furthermore, under the allowed boundary conditions, the Hamiltonian and symplectic 
structure are independent of choice of 12. 

The surface integral §gj. i^w{^, ^2^) will be referred to as the symplectic fiux through 

aE. 

For a diffeomorphism-covariant theory, or a generally-covariant theory on a background 

spacetime, one can show that i/j{$^;(f),5(f)) = 0. Hence the necessary and sufficient condition 
for existence of a Hamiltonian becomes 
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i^uj{^, 62^) = (A30) 

and, furthermore, from the relation (A26) between 50) and B{^; 0), it follows that one 

has 

5(e;0)|aE = M(0))|aE (A31) 

where -6(0) is a locally constructed (n — l)-form. Then on solutions $ the Hamiltonian 
takes the following form: in the case of a diffeomorphism-covariant theory, 

HBi^;^)= i g(e;$)-^^^($) (A32) 

JdT, 

which is a surface integral; and in the case of a generally-covariant theory, 

H{i- $) + $) = / *i^T{^) + / r(e; $) - kB{^) (A33) 

where H{^; $) = Jj. *i^T'($) is the canonical energy associated to $ on E, and Hb{^; $) is 
the surface integral term. 

To conclude, some further features of the Noether charge Hamiltonian will now be de- 
veloped. 

Definition A. 7. The coefficient of an arbitrary compact support variation S(f)\s in the 
equation 

1^ Qs(0, S<l>, C^ct^) - H'^ii- 0, (50) = £H{i] m = (A34) 
yields the Hamiltonian field equations for (j), Sh{C', 0) — 0. 

Theorem A. 8. The Hamiltonian field equations SH{^',(f>) — are equivalent to the La- 
grangian field equations — 0. 

Proof : 

The Hamiltonian (A29) satisfies the variational identity 

Oe(0,50,A0) -^E(e;0,50) = /^^€(^(0)50) (A35) 
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derived from Eqs. (A23) and (A26). Hence, for arbitrary compact support variations S(j)\s, 
£{(t)) = holds if and only if <p satisfies EhU'-, 0) = 0. □ 

A field variation, denoted by 50^/", is a symplectic degeneracy direction if flj:{4>, 54>, 54>j^) — 
holds for arbitrary compact support variations 5(f). Such degeneracies arise whenever the 
Lagrangian 1/(0) admits a gauge symmetry (i.e. a symmetry d^cf) that is locally constructed 
from 0, partial derivatives 9^0, and that depends linearly on a set of parameters x freely 
specifiable as functions on M.) Note that the set {5(1)^} of all degeneracy directions is 
a vector space. Then, a nondegeneracy direction, denoted by is represented as an 
equivalence class in the vector space of all field variations {Scj)} quotiented by all symplectic 
degeneracy directions {5(t)j^}, namely 5(t)D — 5(f)/5(f)j^. This decomposition yields a break up 
of the Hamiltonian field equations (A34) into non-dynamical constraint equations, 

^s(e;0,<^0AA) = O, (A36) 

and dynamical evolution equations, 

H'^{i- 0, 5(l>v) = ^1e(0, 501,, (A37) 

through arbitrary variations ^^a/") ^(t>v with compact support on E. 

Since it assumed that the set {0} of all fields has a hnear (vector bundle) structure, 
the symplectic degeneracy directions 5(f)j^j- can be identified with a corresponding set of field 
components, denoted (j)j\f, which will be called non-dynamical with respect to E. Similarly, 
the nondegeneracy directions 50x> determine a set of equivalence classes of field components, 
denoted 0x>, which will be called dynamical with respect to E. (Note these components 0x) 
and are locally constructed from 0, partial derivatives 9*^0, and spacetime quantities 
associated to E.) Then, from the Hamiltonian variational identity (A35), one can view the 
constraint equations (A36) and evolution equations (A3 7) as arising equivalently through 
the action principle (Al) by variations with respect to 0^ and (f)xi- 

In summary, the Noether charge formalism presented here gives a covariant Hamiltonian 
formulation for Lagrangian field theories in the situation where the underlying time fiow is 
given by a symmetry of the Lagrangian. 
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